Non-Gaussian Statistics of the Soliton Timing Jitter due to Amplifier Noise 
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Abstract 

Based on the first-order perturbation theory of sohton, the Gordon-Haus timing jitter induced by 
amphfier noise is found to be non-Gaussian distributed. Compared with Gaussian distribution given 
by the hnearized perturbation theory, both frequency and timing jitter have larger tail probability. 
The timing jitter has a larger discrepancy to Gaussian distribution than that of frequency jitter. 
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The Gordon-Haus timing jitter in fiber soliton due 
to amplifier noise is usually assumed to be Gaussian 
distributediiSii when the first-order perturbation the- 
ory of soliton^i^ is used. Previous works showed that 
the non-Gaussian timing jitter is induced by soliton 
interactionsLSi^ and regenerationiSiii but not by ampli- 
fier noise alone. When the first-order perturbation of 
soliton is linearized^'^, the Gordon-Haus timing jitter is 
indeed Gaussian distributed. However, if the ordinary 
differential equations from the first-order perturbation is 
not linearized, as shown later, the amplitude, frequency 
and timing jitters are all non-Gaussian distributed. 

In this letter, the amplitude jitter is found to be non- 
central chi-square (x^) distributed, confirming the previ- 
ous simulation results^ and the non-Gaussian statisticaii. 
Both the frequency and timing jitters are also non- 
Gaussian distributed. 

From the first-order perturbation theory of soliton^i^, 
with amplifier noise alone, the soliton parameters are 
evolved according to the following ordinary (or stochas- 
tic) differential equations: 



dA 

dn 

dT 



(1) 
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Not used here, the phase perturbation is not shown in 
(P)-©. All noise terms of n^(^), no (C), and tit (C) are in- 
dependent Gaussian processes with autocorrelation of^-^ 



E{nA{Ci)nA{C2)} 
E{nn{Ci)nn{C2)} 



= Aal6{C, - C2), 



A 



<ylS{Ci ~ C2), 



(4) 
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where a'^ is the noise variance and E{-} denotes ex- 
pectation. From Q to 0, the variances depend on 
the amplitude of A{C,) and amplitude jitter transfers to 
both frequency and timing jitter. If the amplitude in 
the variances of to ® is assumed to be constant 
[A{() fa A{0) = Aq] as a first-order approximation, am- 
plitude, frequency and timing jitters are indeed Gaussian 
distributed. 

The inclusion of amplitude jitter in I^J to © is still 
within the equations of the first-order soliton perturba- 
tion theor y^i^i^ . The nonlinear first-order perturbation 
can be interpreted the repeated usage of the linearized 
first-order perturbation'^. 

Based on and Q), the stochastic differential equa- 
tion (SDE) of the amplitude jitter is dA — ^/AdwA 
where wa is a Wiener process with autocorrelation of 
E{wa{Ci)wa{C2)} = o"^ iiihi(Ci, ^2)- With an initial value 
of A{0) — Ao, the amplitude jitter is 



A, 



1/2 



wa 
2 



(7) 



as a noncentral x^-distributed random process^ii^ii^ with 
variance parameter of jcr^C- The amplitude jitter of Q 
is the solution of a Stranovich but not Ito SDEii^. 
Using (O and ||SJ), the frequency jitter is 



A, 



1/2 



WAiCl) 



dwc 



(8) 



E{nT{Ci)n'T{C2)} = 



12^3 



alSiC, ~ C2), (6) 



where wq is a Wiener process with autocorrelation of 
E{wn{(i)wQ{C2)} = |o'^inin(Ci,C2) and independent of 
the Wiener process of wa ■ 

The timing jitter equation of Q has two terms, the 
first-term of — S7 gives the Gordon-Haus timing jitter in- 
creasing with ajifj second term of 71t(C) is just the 
projection of amplifier noise into timing jitter-'^. The first 
term of Gordon-Haus timing jitter is far more interesting 
than the second term. The SDE of dTcH = —ild^ has a 
solution of 
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Tgh(C) 



(C-Ci)(4/' + 



wa(Ci 



dW: 



(9) 
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Similar to option pricing with stochastic volatihtj*i^, 
the characteristic functions of the frequency ^{C) and 
Gordon-Haus timing jitter Tgh(C) ^re 



(10) 
(11) 



and 



Gi(A) 
G2{X) 



-Ax;(c-ci)^(Ar+^)'rfc 



eu 



(12) 
' } (13) 



where Gi{—X) and G2(— A) are the moment gener- 
ating functions of (^^o^^ + 5'^-4) "^Ci and /o^(C — 

Ci)^ (^A^^ + ^wa) dCi : respectively. 

Based on the Cameron-Martin integrali&, we get 



Gi(A) = 



exp 


-Ma VAtanh 


(% 


-Va) 




cosh^/2 







(14) 



and 

G2(A) 
where 

and 



/a(C) 
/a(0) 



1/2 



exp(AV^Ao/3^-AyloC73), (15) 



1 



/a(Ci) Jo 



(C - C2)VA(C2)dC2 



dCi, (16) 



/a(Ci) = VC^/- 



\/Acr„ 



(C-Ci)^ 



(17) 



where /i/(-) is the iy*'^-oider modified Bessel function of 
the first kind. The function of Gi(A) 114|l is similar to 
the characteristic function of nonlinear phase noiseiii^. 
Using x'^^^ Iu{ax)dx = /^+i(a)/a^^, we get 

(18) 

and 
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Fig. 1. The probability density functions of (a) frequency 
and (b) timing jitter as compared with Gaussian distri- 
bution. 



(3^ = V2(VAa„)-^/2 [^'^"^ '\^'Hl{x)r_\{x)dx. 

Jo * ^ 

(19) 

If the amplitude jitter is approximated as A[C^) « Aq 
in both variances of l@J and we get 



Gi(A)«exp(-AAoC),G2(A) 



exp [ 



(20) 



that are valid for Aq ^ a^C,. With the approximation 
of pO|) . the characteristic functions of H10() and (|ll|l are 
zero-mean Gaussian characteristic functions with vari- 
ance of 



4(C) = ^^oa^C and a|(C) = ^^o^.C', 



(21) 



respectively. Note that the timing jitter has a variance 
increase with C^i. 

From the frequency jitter of ||SJ|, the non-Gaussian dis- 
tribution is induced by the term of ^ WAdwn, i.e., the 
noise and noise interaction. The second-order soliton per- 
turbation also includes noise and noise interaction22i2i. 
However, the equations of to (PI with noise variances 
of to © are directly from the first-order perturbation 
of soliton^. Similarly, the non-Gaussian timing jitter 

of is induced by the term of \ /g^(C — QiiWAdwn that 
also includes noise-noise interaction. 

The probability density functions of frequency and 
Gordon-Haus timing jitter are the inverse Fourier trans- 
forms of the corresponding characteristic functions of 
(|10|l and l|ll|) . respectively. Figs, ^plot the probability 
density of frequency and Gordon-Haus timing jitter as 
compared with the Gaussian distribution. The probabil- 
ity density functions are shown for signal-to-noise ratio 
(SNR) of Afi/{(TlC) = 10 and 20. The horizontal axis 
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Fig. 2. The tail probabilities of (a) frequency and (b) tim- 
ing jitter as compared with complementary error func- 
tion. 

is normalized with respect to the standard deviation of 
frequency cro(C) a-nd timing axiO jitter [see (|2U] for 
Figs, n^a) and (b), respectively. Because the character- 
istic functions of (|10|l and are even real functions, 
the probability density functions are also even functions. 
Figs, njust plot for the positive frequency and timing 
jitter. Comparing the frequency jitter of Fig. ^a) and 
the timing jitter of Fig. Q^b), the frequency jitter has a 
distribution more close to the Gaussian distribution than 
that of timing jitter. 

Figs. 121 plot the tail probabilities corresponding to the 
probability density functions of Figs. ^ The tail prob- 
ability is defined as p{x)dx for a probability density 
function of p{x). The tail probability is compared to the 
complementary error function of ^erfc (x/V2) that is the 
tail probability of the Gaussian distribution in Figs. 
Comparing the tail probability of frequency jitter of Fig. 
|5Ja) and the timing jitter of Fig. I3b), the frequency 
jitter has a distribution more close to the Gaussian dis- 
tribution than that of timing jitter. 

From both Figs. ^ and [3 both frequency and tim- 
ing jitters have a larger tail probability than that of the 
Gaussian distribution. The non-Gaussian distribution 
leads to higher error probability than that of Gaussian 
distribution. 

From both Figs, db) and[2Ib), the distribution of 
Gordon-Haus timing jitter has a large different to Gaus- 
sian distribution at the tail. Comparing the generat- 
ing function of (|15() and the Gaussian approximation 
of H2U|I . the non-Gaussian timing jitter distribution is 

given by the factor of [/a(C)//a(0)]'/' exp(A2a2 Ao/?^) 
that is obviously non- Gaussian. Comparing the gen- 
erating function of (|14|l and the Gaussian approxima- 
tion of H2(J|I . the non-Gaussian frequency jitter distribu- 
tion is given by two factors: the third and higher-order 
powers of tanh(2:) = x — x^ /3 + ■ ■ ■ and the factor of 



sechi/^(C(T„VA/2). 

In conclusion, based on the first-order perturbation 
theory of soliton, both frequency and timing jitters are 
found to be non-Gaussian distributed. Amplitude, fre- 
quency, and timing jitters are all Gaussian distributed if 
the equations from perturbation are linearized. Without 
linearization, the noise projected into frequency and tim- 
ing jitters are modulated by the amplitude jitter, lead- 
ing to non-Gaussian distribution. The timing jitter has 
larger discrepancy to Gaussian distribution than that of 
frequency jitter. 
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